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SUMMARY 


Methods of controlling the trajectories of high-drag — low-lift 
vehicles entering the earth’s atmosphere at angles of attack near 90° 
and at initial entry angles up to 3° Q:ce studied. The trajectories 
are calculated for vehicles whose angle of attack can be held con- 
stant at some specified value or can be perfectly controlled as a func- 
tion of some measured quantity along the trajectory. The results might 
be applied in the design of automatic control systems or in the design 
of instruments which will give the human pilot sufficient information to 
control his trajectory properly during an atmospheric entry. 

Trajectory data are compared on the basis of the deceleration, 
range, angle of attack, and, in some cases, the rate of descent. The 
aerodynamic heat-transfer rate and skin temperature of a vehicle with 
a simple heat-sink type of structure are calculated for trajectories 
made with several types of control functions. 

For the range of entry angles considered, it is found that the angle 
of attack can be controlled to restrict the deceleration down to an 
arbitrarily chosen level of 3g. All the control functions tried are 
successful in reducing the maximum deceleration to the desired level. 
However, in order to avoid a tendency for the deceleration to reach an 
initial peak, decrease, and then reach a second peak, some anticipation 
is required in the control function so that the change in angle of 
attack will lead the change in deceleration. When the angle of attack is 
controlled in the aforementioned manner, the maximum rate of aerodynamic 
heat transfer to the skin is reduced, the maximum skin temperature of the 
vehicle is virtually unaffected, and the total heat absorbed is slightly 
increased. The Increase in total heat can be minimized, however, by 
maintaining the maximum desired deceleration for as much of the trajec- 
tory as possible. From an initial angle of attack of 90®, the angle-of- 
attack requirements necessary to maintain constant values of deceleration 
(lg to 4g) and constant values of rate of descent (450 to 1,150 ft/sec) 
as long as it is aerodynamically practical are calculated and are found 
to be moderate in both magnitude and rate. Entry trajectories made with 
these types of control are presented and discussed. 

*Title, Unclassified. 
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INTRODUCTION 


successful recovery of satellite and near-satellite vehicles 
it is generally necessary to maintain deceleration and heating within 
certain limiting values. In addition, it is desirable to be able to 
predict with reasonable accuracy the trajectory and range of the vehicle. 
A logical means of achieving these objectives is to control the trajec- 
tory by varying the aerodynamic forces generated during entry into the 
atmosphere . 

At present, the National Aeronautics and Space Administration and 
other research agencies are investigating factors important to the devel- 
opment of satellite vehicles as well as a number of general satellite con 
figurations which hold promise of achieving desirable and predictable con 
ditions during their return to earth. This paper deals with the control 
of the point-mass trajectories through the earth's atmosphere of one 
general class of space vehicles, namely, vehicles which have a large 
flat-plate area which is maintained nearly normal to its flight path 
during most of the entry. This class of vehicles includes circular- or 
delta-wing gliders as well as capsules having blunted faces. Methods 
'-cntrolling the deceleration and trajectory include either varying or 
maintaining constant the angle of attack in the range near 90° . By 
maintaining the angle of attack near 90°, a strong detached shock wave 
is produced ahead of the vehicle which in turn dissipates a large per- 
centage of the energy into the atmosphere with a consequent lower heat 
transfer to the vehicle itself. By changing or maintaining the angle 
of attack to values smaller than 90°, a component of the resultant 
aerodynamic force is produced normal to the flight path. This lifting 
force reduces the rate of descent into the denser atmosphere and conse- 
quently the deceleration and heating that would result in this denser 
atmosphere . 

The first consideration, that of maintaining a large ratio of wave 
drag to friction drag, is an important result reported by Allen and 
Eggers in reference 1. The second consideration which involves the 
production of small lift-drag ratios to control point-mass trajectories 
has been widely recognized as a most logical approach to trajectory 
control . Some of the effects of a small constant lift-drag ratio during 
atmospheric entries are given 'in reference 2. The purpose of this paper 
is to utilize these concepts and to determine the effects of controlling 
the angle of attack, and hence the lift, as a function of the measurable 
trajectory variables in such a manner as to achieve low levels of 
deceleration which will not incapacitate a human pilot and to minimize 
the increase in the total temperature and the range compared with that 
of the zero-lift case for the same vehicle. The results might be 
applied in the design of automatic control systems or in the design of 
instruments which would give the human pilot sufficient information to 
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control his trajectory properly during the entry. As an aid in 
achieving these ends, studies are also made of families of trajectories 
in which the angle of attack is varied to maintain one of the flight 
variables constant, namely, either deceleration or rate of descent. 

Although the assumptions used in the trajectory calculations are 
believed to be adequate to provide quantitative results, the calcula- 
tions of temperatures are only qualitative since the temperature over 
the surface of the vehicle will depend on many specialized factors of 
design and material. 


SYMBOLS 


a- 


n 


a v 




D 

F 


H 

k 

K-i 


deceleration in direction of resultant aerodynamic force, -a z 

component of deceleration along flight path, V/g 

deceleration along Z-axis of vehicle 

speed of sound, (dp/d 0 )/ 

' ** H isentropic 

heat capacity of material 
aerodynamic force coefficient, 2F/pV^S 

drag, F sin a 

resultant aerodynamic force 

acceleration due to gravity 

radial height above surface of earth, r - r e 

constant for exponential variation of density 
emissivity -absorptivity factor 

variation of angle of attack with deceleration (see eq. (20)) 

variation of angle of attack with rate of change in decelera- 
tion (see eq. 21)) 


variation of angle of attack with temperature 
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L 

m 

M 

P 

<1 

Q/A 

r 

r* 


R 

As 


S 


t 

T 


T 1 


Tj* 


T 


2 


V 


V 


t 


W 

w/s 

Ui/A 

X 


lift, F cos a 
mass of vehicle, W/g 
Mach number, V/c 
free-stream pressure 
free-stream dynamic pressure, 

heat-flow rate per unit area 

radial distance measured from geographical center of earth 

effective radi us of total surface area of vehicle normal to 
flow, J S/jt 

perfect gas constant for air 

perimetric distance over surface of earth traversed since 
zero time (range) 

frontal area of vehicle (assumed to be equivalent to wing area) 
time 

free-stream temperature 

vehicle skin temperature based on WpC/A of 0.8 
vehicle skin temperature based on WpC/A of 0.2 
radiation equilibrium temperature 
free-stream velocity 

tangential velocity, component of total velocity normal to 
radial vector 

weight of vehicle 

wing loading of vehicle 

weight per unit surface area of vehicle skin material 
vehicle body axis alined with wing-chord plane 
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Z vehicle body axis normal to wing -chord plane 

a angle of attack 

P(t) an explicit function of time used in appendix A 

7 flight-path angle 

7 0 initial flight-path angle , entry angle 

ratio of specific heats for air 

5 angle between velocity vector and reference horizon 

© angle between reference horizon and tangential velocity, 

angular rotation of radial vector to point mass 

^ 1^2 temperature constants 

p free-stream density 

P e constant for exponential variation of density 

a Stef an -Boltzmann constant for radiation 

Symbols used only in appendix C: 
aQ horizontal component of deceleration 

u horizontal component of velocity 

u = u jj gr 2 v/^gr 

Z nondimens ional parameter of reference 2 

0 — l/H 

9 = 7 

M) = () n -( 

Subscripts : 

o initial condition, that is, value at zero time 

e referenced to mean surface of earth 

d design 
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Derivatives with respect to time are denoted with dots over the variable. 
For example, dy/dt = 7. 

EQUATIONS OF MOTION 


This study is concerned with entry into the earth's atmosphere of a 
vehicle which is traveling at or near the velocity of 26,000 feet per 
second and is entering the atmosphere on a flight path which initially 
differs from the instantaneous horizon by a small angle (-7) not exceeding 
3°* On the basis of current knowledge of the atmosphere the trajectory 
studies were started at an initial altitude of 350,000 feet. It is 
assumed that above this altitude the effect of the atmosphere on the 
flight path is negligible.^ A schematic diagram of a typical trajectory 
showing this initial condition is shown in figure 1. The forces, acceler- 
ation, velocities, and displacements are shown in their positive sense 
except as noted. 

The vehicle considered in this study is assumed to have a large 
flat-plate area such as a wing which is stabilized nearly normal to the 
flight path (that is, at an angle of attack of nearly 90°) at the time 
°f the initial conditions . It is further assumed that this angle of 
attack may be maintained at a constant value or changed without lag to 
any desired value during its trajectory within the atmosphere. The 
resulting motions are assumed to be those of a point mass and no moment 
equation is included in the equations of motion. 

The governing equations of motion, when the notation of figure 1 
is used, are as follows: 


mV = -D - W sin 
mVS = L - W cos 7 
h = V sin 7 
A V 

0 = — cos 7 
r 


( 1 ) 


For the relatively small values of altitude covered, the accelera- 
tion due to gravity and the distance from the vehicle to the center of 
the earth are assumed to be constant; that is, 

1 The effect of the atmosphere on vehicle temperatures at these alti- 
tudes however is not negligible. (See the section entitled "Temperature 
Equilibrium Equations . " ) 
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g = 32.2 ft/sec 2 
gr = 25,863 ft/sec 

An exponential variation of density with altitude is used such 

that 

p(h) = P e e' h / H (2) 

with 

p e = 0.003 
H = 23,000 ft 


Newtonian flow is assumed at the surface of the flat-plate area 
during the entire trajectory although this assumption will undoubtedly 
be violated during the final stages of the trajectory where the 
velocity falls below 5 >000 feet per second. Accordingly, a force 


F = |pV 2 SC F = -mga z (3) 

is assumed to act normal to the flat-plate area at all times. The com- 
ponents of this force along and normal to the flight path are denoted 
as the drag 


D(a) = F sin a 


oo 


and the lift 


L(a) = F cos a 


(5) 


in the usual sense. A value for the aerodynamic-force coefficient normal 
to the plate was chosen as 


Cp = 1.7 
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and was assumed to be invariant with angle of attack over the range of 
conditions covered. This value of C;p is slightly less than the pure 
Newtonian pressure coefficient of 2.0 and slightly greater than the 
(asymptotic) value of about 1.6 which is predicted for the pressure 
coefficient on a cone of maximum apex angle for an attached shock at 
hypersonic Mach numbers. (See ref. 3 -) 

The variable 5 is related to © and y by 


5 = 7 -© ( 6 ) 


as defined by the sign convention of figure 1 . By substituting equa- 
tions (2) to (6) into equations (l), the equations of motion may be 
written as follows: 


— - -a^ sin a - sin 7 


V7 / y 2 \ i 

— = cos a + - lycos 7> 


h = V sin 7 


x n 


-a^ = 


'p£A-h/H v 2 

2 W/S 


( 7 ) 


( 8 ) 


where is the acceleration measured normal to the flat plate, as is 

shown in figure 1. 

The range on the surface of the earth is given by the integral 

f 0 d) r t 

As = r e / d© = r e / cos 7 dt (9) 

Jo Jo 

Since r was assumed to be constant and equal to r 0 , no distinction is 
made between the range on the surface of the earth and that at altitude h. 

When the angle of attack and the three initial conditions are speci- 
fied, equations (7) to (9) define the trajectory from these initial 
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conditions. Another type of solution may be found when some variable 
or combination of variables other than angle of attack is specified and 
the equations are solved for the variation in angle of attack necessary 
to produce this condition. Both types of solutions are investigated in 
subsequent sections. 


TEMFERATURE EQUILIBRIUM EQUATIONS 


For the trajectories where temperature calculations are presented, 
a simple model was used to represent a unit area of the exposed surface 
of the vehicle . The rate of heat flow into the material from the boundary 
layer was calculated by using an expression derived in reference A for 
the stagnation point of a rounded blunt body. This expression, given in 
units of Btu/sec-ft 2 , is 

(;) „ . - oWW-5 do) 

\ ' aerodynamic \ / 


where the constant 0.010 is determined by the roundness of the surface, 
r* is the effective radius of the total surface area normal to the flow, 
and p and M are, respectively, the free-stream static pressure and 
Mach number. The constant 0.010 is somewhat less than the value used 
by Romig since it is assumed that the curvature of the undersurface of 
the vehicle is less rounded than that of the model used by Romig. The 
value of r* was taken as 8 feet which corresponds to a vehicle flat- 
plate area of 200 square feet. 

The net rate of heat transfer by radiation from the surface at a 
temperature T-^ to the surroundings at temperature T2 was calculated 

from the expression: 


r) = ka 

^ 'radiation 


{?l k - T 2 k ) 


( 11 ) 


where a is the Stefan -Boltzmann constant, 0.48 X 10“12 — — and 

ft 2 -sec-°R 4 

k is the emissivity -absorptivi ty factor. Unpublished calculations and 
measurements indicate that at an altitude of about 100 miles the radia- 
tion equilibrium temperature between the surface and the surroundings 
will be about 58O 0 R. Since the calculated aerodynamic heating rate 
(based on the assumed model atmosphere) is negligible at this altitude, 
the initial equilibrium temperature of the vehicle will also be about 
580° R. 
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It was assumed that the surface radiated from only one side but at 
an effectiveness of unity (k = 1 . 0 ). The result of the two assumptions 
is more realistic than the results with either assumption singly. 

The heat absorbed by the structure is dependent on the product of 
the weight per unit area and the heat capacity of the material used. 

For the calculations herein, the rate of heat absorbed per unit area 

W-.C . 

r T, hi 


was calculated for l/8-inch solid beryllium such that 

W-, 

-£■ = 1.1844 lb/sq ft 


c = 0.675 P tu 

lb-°R 


W-i c 

The product is thus 0 , 8 , In some cases^ the material was con- 

sidered to have less heat capacity (for example, a thinner gage material) 
and the product of the two terms was reduced to 0.2. However, unless 
otherwise specified, the value 0.8 was used in all temperature 
calculations . 

A third element of heat transport, conduction between the skin 
material and the internal structure, is neglected in these calculations. 
One further simplification is employed by assuming that this unit area 
of material is always located at the stagnation point (or region) of the 
vehicle's surface in spite of the fact that the angle of attack is being 
varied. Some justification for this assumption may be found in 
reference 5. 

The heat capacity of the material together with the two elements of 
heat transport considered, namely, conduction from the boundary layer to 
the material and radiation between the material and its surroundings 
constitute a balance within the chosen control surface. Thus, 


^ — ■ ■ , - 

For the purposes of this paper, the material need not be an integral 
part of the structure but may be a temperature sensor which is insulated 
from the rest of the structure. 
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a'j - (a) = (13) 

A /aerodynamic ' A 'radiation A 


throughout the trajectory. When the indicated substitutions and grouping 
of the coefficients are made, the governing differential equation becomes 


T l + + l 2 M5 ' lp0 ' 5 


where 

t _ koA _ q. 6 x io -12 

1 cw L 


= 0- 01A = 0.004417 

CW 1 (r*) 0,5 


Unless indicated otherwise, equation (l4) was used with these values of 
S -|_ and | 2 • 

At each position on the trajectory the pressure and Mach number, 
were calculated from the altitude and velocity data from the expressions 

p = Rp(h)T(h) ( x 5) 

M = I = V— (16) 

c y 7l RT(hy 


where 


and R were treated as constants having the values 


R = 1716.6 ft 2 / sec-°R 
7]_ = 1*^ 


The variables p and V were taken directly from the trajectory data, 
whereas the free-stream temperature T was computed by using straight- 
line segments to approximate the altitude -temperature plots of reference b 
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At 350,000 feet (h Q for most of the trajectory studies) the 
equilibrium surface temperature calculated by equation (l4) exceeds 

1.000 R and the assumption that Tq is zero at that condition is not 

justified. Therefore, trajectories which include temperature calcula- 
tions were initiated at 550,000 feet where aerodynamic heating is negli- 
gible, Tq ~ 0, and the surface temperature of the vehicle and the 

radiation equilibrium temperature are about equal (580° R). Since the 
angle of attack has virtually no effect on the trajectory above 

350.000 feet, only a few cases were calculated at 550,000 feet to 
establish the value of Tq and Tq at 350,000 feet for each initial 

condition needed. 


TRAJECTORY STUDIED 


In order to provide insight into means of regulating the maximum 
deceleration during the entry, the conditions governing the maximum 
deceleration are derived. The necessary condition that the deceleration, 
given by equation (8) be constant or at least an instantaneous maximum 
is imposed by setting its time derivative equal to zero. Thus for 
these extremum, the condition 

= 57i e ' h/H (-l v2 - 2rt ) = 0 


leads to a necessary condition between the rate of descent and the time 
rate of change of velocity; that is. 


V = Vh_ 
S 2gH 


(17) 


Expressions for these two derivatives everywhere on the trajectory are 
given by equation (7). Substituting these expressions into equation (17) 
and rearranging the terms yield the desired expression for the maximum 
value of the deceleration 


^ an ^max 


sin 7 / 
sin a\2gH 



(18) 
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For the trajectories controlled with aerodynamic lift, equation (l 8 ) 
shows that the maximum value (or some desired constant value) of deceler- 
ation depends directly on the instantaneous values of velocity, angle 
of attack, and flight-path angle (or rate of descent since h = V sin 7 ). 
As long as the angle of attack is near 90°, its variation at the instant 
of consideration has very little effect on the magnitude of the decelera- 
tion. However, if the angle of attack is used to produce a force normal 
to the flight path for a sufficient time prior to the time of the maximum 
deceleration, the magnitudes of both flight-path angle and velocity at 
the time of maximum deceleration can be reduced to some degree. Likewise 
if this same force is controlled with the angle of attack throughout the 
entire atmospheric portion of the entry, it may be aerodynamically 
possible to hold the deceleration, heating, and range to within some 
prescribed limits. The angle -of -attack requirements and some methods of 
application for the aerodynamic control of entry vehicles are investigated 
in the following sections. 

The trajectory results are obtained for a vehicle having an assumed 
value of the ratio of wing loading to force coefficient W/CpS; however, 
it is shown in reference 2 (and in appendixes A and C of tnis paper) that, 
with an exponential altitude -density relationship, the results are appli- 
cable to any value of W/CpS if the altitude corresponding to each point 
on the trajectory is shifted by an amount proportional to the logarithm 
of W/C F S. 


Constant-Deceleration Trajectories 

The equations of motion given previously may be used to solve for 
the variation of angle of attack and the resulting variations of altitude 
and velocity necessary to maintain a constant deceleration during a por- 
tion of an entry into the earth’s atmosphere. The entry trajectory starts 
from a zero-g condition and the deceleration increases as aerodynamic drag 
is encountered. During the final stages of the entry, the deceleration 
again may decrease as the dynamic pressure falls to values too low to 
sustain a prescribed value of deceleration. The portion of the trajec- 
tory involving constant deceleration is therefore accompanied by 
transition periods of lower deceleration at the start and finish of the 
entry. The constant-deceleration period is of interest, however, because 
a constant-deceleration trajectory will allow the shortest path through 
the atmosphere without exceeding deceleration values imposed by structural 
or human physical limits. Furthermore, the oscillations or skipping 
motions which accompany some constant L/D trajectories (see ref. 2) are 
eliminated. The results may also prove to be useful as an ideal guide 
for other types of control. 
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Since the maximum deceleration occurs at all times colinear with 
the resultant force , the calculation is restricted to the aerodynamic 
force and, hence, to the deceleration or -a z . In order to main- 

tain a constant value of deceleration for a given vehicle, the dynamic 
pressure must remain constant at its corresponding value and, when the 
exponential approximation to the atmospheric density is used, the 
relationship between altitude and velocity is fixed. When equation (8) 
is rearranged, this relationship becomes 


h = H log e 


C F p e v£ 
2W/S ^ 


( 19 ) 


A plot of the variation of h with V for various values of is 

shown as figure 2 for a vehicle with a wing loading of 20. It may be 
seen that from the initial high-velocity high-altitude condition a 
large reduction in velocity would occur with only a small reduction 
in altitude. At velocities below 5; 000 feet per second a large reduc- 
tion in altitude would occur with relatively small changes in velocity. 

In order to illustrate some of the differences between two types 
of trajectories, data from a typical high-drag nonlifting (a = 90°) 
entry trajectory are superimposed on the lines of constant deceleration 
in figure 2. For the zero lift case the flight-path angle is -3° at 
350,000 feet and time is indicated on the figure by using dots at 
10-second intervals. The vehicle enters the atmosphere at a velocity 
which has changed very little by the time the deceleration reaches 1 g. 
Thereafter, the vehicle follows a high deceleration path which dissipates 
most of the energy over a relatively small range of altitudes. 

The initial portion of the nonlifting trajectory is of particular 
interest. Above 320,000 feet the dynamic pressure is less than 1 pound 
per square foot. If the descent is not checked, a 3g deceleration is 
reached at 230,000 feet, 56 at 220,000 feet, and so on. This means that, 
during the time required to descend 80,000 to 100,000 feet, control of 
the trajectory by means of the aerodynamic force is available and 
necessary if the vehicle is to -be guided toward an asymptotic inter- 
section with a constant 3g to 56 deceleration trajectory. 

The conditions necessary for a vehicle to remain on a constant- 
deceleration trajectory are investigated next. Development of a form 
of the basic equations suitable for calculating constant-deceleration 
trajectories is given in appendix A. The equations were used to compute 
the trajectories and the variations in angle of attack necessary to achieve 
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the trajectories for the cases of = Ig, 2g, 3g ; and 4g . The computa- 


tions were performed on a digital computer with a triangular (two-point) 
integration process. 5 Each trajectory was started at the flight condi- 
tions V, 7, h which satisfied the equations when a was exactly 
90°, that is, the maximum drag condition, (See eqs . (Al6) and (A17).) 


Time histories of the trajectories are presented in figure 3. It 
is noted that the initial value of h is almost linearly dependent on 
the magnitude of the deceleration. From its initial value of 90° the 
angle of attack decreases; this decrease in angle of attack Increases the 
lift, prevents a rapid increase in flight-path angle or rate of descent, 
and thus maintains a constant deceleration. As the velocity decreases, 
however, the requirement (eq. (19)) for lower altitudes leads to ever- 
increasing rates of descent and flight-path angles. Except in the 1 g 
case, this demand for higher rates of descent eventually leads to angles 
of attack greater than 90° and, hence, a negative lift is required to push 
the vehicle down into the denser atmosphere if a constant deceleration is 
to be maintained. 


The time histories were discontinued at velocities less than 
4,000 feet per second. It may be expected, however, that a and -7 
continue to increase and the velocity continues to decrease until the 
desired value of deceleration can no longer be obtained. The exception 
is the 1 g trajectory which the equations indicate can be maintained to 
sea level where the terminal velocity is about 100 feet per second for 
a = 90° and w/S = 20. Since the assumed aerodynamic characteristics 
become inaccurate at low angles of attack as well as at low velocities, 
the data given for the 1 g trajectory should be viewed as approximate. 

From the trajectory data of figure 3j the variation of angle of 
attack and flight -path angle with velocity is crossplotted in figure 4 
to eliminate time as the independent parameter. The results, presented 
in this form, are very general since wing loading, time, altitude, and 
force coefficient are not involved in these data as may be seen from the 
equations of appendix A. The effect of the aforementioned parameters is 
contained entirely in the altitude -velocity relationship given by 
equation (19). Therefore, if W/S or Cp are changed, the variables 

7, a, h, tf, and all of their time derivatives remain unchanged as 
functions of velocity but each velocity on the trajectory occurs at a 
different time and altitude. 


^By the use of certain approximations, analytical solutions which 
are valid for a large part of each trajectory may be obtained. These 
closed-form solutions and the corresponding assumptions are given in the 
latter part of appendix A and in appendix C. 
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Control of Angle of Attack As a Function of Deceleration 

In order to reduce the peak deceleration from the value obtained 
with a nonlifting entry trajectory, the control of angle of attack as a 
linear function of some trajectory variables was tried. One obvious 
quantity to use for this purpose is the deceleration measured along the 
Z body axis of the vehicle. Since the resultant aerodynamic force in 
Newtonian flow Is also along this axis, no resolution of accelerations 
is necessary. For this type of control 


a = a 0 - 


( 20 ) 


and a no-lag response in a to a measured value of is assumed. 

This would be an Ideal choice if, by adjusting K]_, the deceleration 

could be made to reach some predetermined maximum value and remain in 
close proximity to that value until the low-speed (that is, under 
5,000 feet per second) condition is reached. The constant-deceleration 
data of figure k could then be used as a guide in the choice of K^. 

In order to illustrate these points, the time history of an entry 
trajectory with the angle of attack held constant at 90° is presented 
in figure 5(a) and an entry with a = 90° - 3a n is shown in figure 5(^0 • 

The trajectory equations are equations (7), (8 ), and (9) and the initial 
conditions in both cases are 


a o = 90° 

To " - l0 

h Q = 550,000 ft 
V 0 = 25,863 ft/sec 


When the angle of attack Is held constant at 90° (fig. 5 ( a ))j the 
deceleration rapidly builds up to a peak value of 8g. With the angle of 
attack varied according to equation (20), figure 5(b) shows that the rate 
of descent is reduced and the maximum deceleration is less than kg. 

In the constant-deceleration data of figure k the minimum angle of 
attack for a constant 3g entry is 78° and this occurs at a velocity near 
11,000 feet per second. A gain K-j_ = k in equation (20) would likewise 
produce an angle of attack of 78° whenever Therefore, the 
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relative effect of increasing the gain to 4 is illustrated in figure 6. 
Time histories of deceleration and angle of attack are compared in 
figure 6 for values of Kp of 0 , 3 > and 4 . It may be seen that the 
deceleration was held in the vicinity of 3g when the gain was 4 although 
some overshoot did occur. The velocity of 11,000 feet per second 
occurred at 365 seconds for the Kp = 4 case and occurred at 345 seconds 

for the Kp = 3 case . 

In figure 6 and in most ensuing figures, the range As is given for 
each trajectory. The range is defined as the distance in miles traveled 
over the surface of the earth (nonrotating) from the initial position in 
space to a point on the trajectory where the range is increasing at less 
than 100 feet per second. At this condition total velocity is under 
5,000 feet per second and the flight path is very nearly vertical 
(7-^.- 90°). Since for all trajectories the variations of As, V, h, 
and 7 are nearly identical to those given in figure 5> they are not 
included in the figures . 

Figure 7 shows the effect of initial flight-path angle on the 
deceleration and rate of descent when the angle of attack is controlled 
according to equation (20). When the initial flight-path angle is 
between 0° and -1°, the deceleration reaches about the same maximum value. 
However, at entry angles greater than 1°, the increased rate of descent 
produces' a more rapid increase in deceleration and the consequential 
reduction in angle of attack does not have time to reduce the flight - 
path angle and rate of descent sufficiently before large decelerations 
have been reached. Thus it may be seen from figure 7 that this type of 
controller, which may be efficient at some specified entry angle, does 
not possess the ability to anticipate the more rapid changes in 
deceleration which occur at larger entry angles or higher rates of 
descent. The key here, then, is anticipation. Several means of 
obtaining the necessary anticipation are available such as measurements 
of 7, h, a n , and Tp. If it is known during the early stages of the 

entry that 7 is large, the angle of attack can be reduced to some value 
less than 90° before the deceleration starts to build up. (This change 
presumably would require reaction controls . ) The same thing could be 
done by measuring ft in the early stages of the entry. Further consid- 
eration will be given to these and temperature as methods of control in 
subsequent sections of the paper. 

The third method of anticipation, that of a^, can be easily added 

to the present study. The angle of attack and, hence, the lift and 
flight path may be changed earlier in the entry by making 


a 


a o - K ian - K-2&n 


( 21 ) 
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and. by properly choosing K^. As in previous cases, a perfect response 
in angle of attack to the controls is assumed. 

The complete trajectory of a vehicle with this type of angle -of - 
attack control is shown in figure 8 for an initial flight-path angle of 
-1°. Gains were set at Kp = 3 deg/g, and Kg = 250 deg/g/sec based 

on measurements of a n and a n in the preceding cases. In figure 9 
the relative effect of this type of control is compared with the 

control and the constant 90° angle -of -at tack entry of figure 5* The 
effect of the anticipation term in the angle -of -attack control is to 
reduce notably the magnitude of the first peak in the deceleration 
while the magnitude of the second peak remains virtually the same. 

A range of values for K2 is illustrated in figure 10. With the 
initial conditions of a Q = 90°, y Q = -1°, and Kp = 3 common to all 

three cases, the variation of deceleration (to an enlarged scale) and 
angle of attack with time is shown for the cases of K2 = 150* 25O, and 
350. The figure illustrates that increasing the value of K2 primarily 
suppresses the rate of increase in the deceleration. 

Larger values of the initial flight -path angle present a more 
severe test of the angle -of -attack controller and trajectory data at 
several entry angles are presented in figure 11. Gains of Kp = 4 and 

Kg = 250 and an a 0 = 90° are used in the angle -of -at tack control 

while entry is initiated at flight-path angles of -1°, -2°, and -3°. 

The deceleration, rate of change of altitude, and angle of attack of 
figure 11 may be compared with those of figure 7 ^ where no anticipation 
was used. With anticipation the deceleration was held within the values 
of 3-^g to 2.1g. In order to obtain these lower values of decelerations 
over the range of entry angles, it is necessary to make relatively large 
changes in angle of attack but at a relatively low rate. For example, 
when 7 0 = -3° the minimum angle of attack is 59° "but the maximum rate 
of change of angle of attack required is 1° per second. These values 
compare with the minimum angle of attack of 67.2° and a maximum rate of 
O.52 0 per second obtained in the y Q = -3° case of figure 7* 

The effect of relatively large changes in angle of attack on the 
deceleration, heating rates, and the temperature of the vehicle is shown 
in figure 12. In this figure a third possible control function is used, 

^The angle of attack, although not shown, is easily computed from 
figure 7 since a = 90° - 4a z . At the maximum decelerations, angles of 

attack of 77-2°, 73-8°, and 67.2° are obtained for the respective entry 
angles of -1°, -2 , and -3°. 
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namely, the case where is set to 0 and a = 90° - 250a n . Trajec- 

tories using this control function are compared with trajectories when 
a = 90 ° on the basis of deceleration, angle of attack, aerodynamic 
heating rate to the vehicle, and skin temperature for initial flight- 
path angles of -1°, -2°, and -3°. When compared with the data of the 
preceding figure (fig. ll), the effect of changing K-j_ from 4 to 0 is 

to increase at all three entry angles the value of the second peak in 
the deceleration by about 1.5g. From a comparison of the a = 90° - 250a n 

trajectories with the a = 90 ° trajectories, the following conclusions 
are drawn. When the angle of attack is changed in proportion to a n , the 

time histories of deceleration are very similar and have almost the same 
maximum value for entry angles up through 3°* Heating rates and maximum 
skin temperatures are virtually unaffected by large reductions in angle 
of attack from 90 ° compared with the value obtained when the angle of 
attack is held constant at 90°. Note that the skin temperature and 
heating rate reach their maximum values before the deceleration reaches 
a maximum (exception: a submaximum in the y 0 = - 3 ° case); thus, it 

appears that a reduction in angle of attack to values less than 90 ° 
during the initial phase of the trajectory produces as its principal 
effect a less rapid decrease in the heating rate and temperature during 
the final phase of the trajectory. This effect results in a slightly 
larger total heat which must be absorbed by the vehicle but this increase 
occurs at temperatures less than the maximum. 

In all the preceding trajectories it has been assumed that the 
vehicle was at an angle of attack of 90 ° when entry was initiated (that 
is, by the time the vehicle reached 350,000 feet). In figure 13, with 
an entry angle of - 1 °, the entry is initiated at angles of attack of 70 ° 
and 110° with the a Q = 90° case reproduced for comparison. The type 
of control assumed for this study is one in which 


a = a o - *1*0. - K 2 *n 


where now 


K i 


a o - a d 
a n,d 


In the expression for K-j_, a n ^ I s the "design deceleration" taken as 

3 g and is the corresponding minimum angle of attack ( 78 °) for a 

3 g deceleration as obtained from figure 4. Thus with this control law 
the angle of attack varies from the value a Q at the start of the entry 
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(where the deceleration is zero) to the value of ct,^ when the deceler- 
ation reaches design value When a Q is greater than a<j, 

is positive and, as shown in previous figures, a positive produces 

a stable variation of angle of attack with deceleration for entry trajec- 
tories. When a Q is less than a^, is negative, a condition which 

has not been investigated previously. 

For the purposes of computation, the values of K L were rounded 
off to the nearest integer; thus, 


a Q = 110° K x = 10 

a Q = 90 ° K L = 4 

<x 0 = 70° K x = -3 

As shown in figure 13, the control was effective even when is 

negative. The effects of changing density together with the stabilizing 
contribution of the term in the controller act to nullify any rapid 

increase in angle of attack which would in turn increase the decelera- 
tion beyond the desired limits. When this type of controller is used, 
the differences in range due to the initial difference in angle of 
attack are very pronounced: With reference to the a Q = 90° case 

the range is decreased by about 250 miles when a Q = 110 ° and increased 
by about 900 miles when a Q = 70 °. 


Trajectories for Vehicles Trimmed at Angles of 
Attack Less Than 90° 

In the preceding sections it has been shown that the deceleration 
during reentry may be limited to small oscillations about some mean value 
when the angle of attack is reduced from 90 ° to some lesser value as a 
function of an and a^. When the angle of entry was greater than 1 °, 
however, large changes in angle of attack were necessary to prevent large 
overshoots in deceleration. By entering the atmosphere at angles of 
attack less than 90 < “ l j it may be possible to obtain a compromise between 
drag and deceleration which would allow smaller and even less rapid 
changes in angle of attack and still retain a high ratio of wave drag to 
viscous drag in order to minimize heating of the vehicle. In figure 14 
trajectory data for four entry angles are shown for the condition where 
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the angle of attack is held constant at 79° • 'For entry angles up to 2° 
the deceleration was limited to 3-3g while at a 3° entry the deceleration 
reached 4.8g. However, these values represent an appreciable reduction 
in deceleration over the 8g to llg range when the angle of attack is held 
constant at 90° during the entry. (See fig. 12.) 

A comparison of the aerodynamic heat-transfer rates and vehicle 
skin temperatures is shown in figure 15 for the a = 90° and a = 79° 
cases and for entry angles of 1° and 3°* When the entry angle was 1°, 
the vehicle maintained at a = 79° had a lower aerodynamic heat-transfer 
rate and a maximum skin (WpC/A = 0.8) temperature that was lower by 

100° R than that for the vehicle maintained at a = 90° • Total heat 
transferred to the vehicle through aerodynamic heating, as measured by 
the time integral of the aerodynamic heating rate, appears to be one- 
third greater when the angle of attack is held at 79°. 

The 7 0 = -2° case is not shown in figure 15 since it afforded 
no additional information. When y 0 - -3° the maximum values of aero- 
dynamic heat-transfer rate and skin temperature are about the same at 
either angle of attack. For both entry angles shown, the total heat 
transferred is considerably increased when a = 79° as compared with the 
a = 90° case. The increase in heat transfer may be further compared 
with the cases shown in figure 12. In figure 12 the angle of attack was 
increased to angles greater than 90° as soon as the deceleration had 
been checked, rather than holding it constant at 79° throughout the 
trajectory as was done in figure 15. It may be seen that the type of 
angle-of -attack control used in figure 12 results in a total heat transfer 
which is much lower than that obtained in figure 15 when a is held con- 
stant at 79° • This difference in total heat transfer may be important for 
vehicles which depend to some degree upon ablation materials for heat 
protection. 

A comparison between entries made at a constant angle of attack 
(other than 90°) and entries made with the angle of attack varied (from 
its initial value) as a function of is shown with two examples in 

figure l6. In figure 16(a) the time histories of and a are shown 

for an initial flight-path angle of 3° and with a = 79° - lOOa^. The 
effect of the anticipation is to reduce the maximum value of by 

about 1.2g with an 11° variation in a and to increase the range by 
104 miles. The rate of change of angle of attack did not exceed O.5 0 
per second. 

In figure l6(b) trajectory data are shown for the same conditions 
as figure l6(a) except that a 0 is 84°. When the angle of attack is 

held at 84°, the deceleration reaches 6.8g. When a is changed in 
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proportion to a,^ a 12° variation in a is sufficient to reduce the 
deceleration by 2.4g while the range is increased by 110 miles. 

From a number of such trajectories a summary plot of range and 
maximum deceleration as a function of entry angle is presented in 
figure 17. It was assumed in obtaining these data that 79 0 was the 
desired angle of attack but that a ±5° error in a D might occur. 

Therefore , the data are given for initial values of a of 84°, 79°, 
and 74° for cases where the a n gain is 0 and 100. The data show that, 
when the a n control is included, the deceleration is reduced, particu- 
larly at the higher entry angles, while the range is extended by less 
than 110 miles. 

Within the conditions investigated, the data of figure 17 indicated 

that 


Range <x 7 0 "®*5® 

when the anticipation control is included. The range does not exactly 
follow a power relationship, however, when the angle of attack is held 
constant. 


Trajectories Made at a Constant Rate of Descent 

.In several places in the preceding sections it has been pointed out 
qualitatively that rate of descent and maximum deceleration are closely 
related during the atmospheric entry. Allen and Eggers (ref. l) have 
shown that, for a nonlifting body passing through the atmosphere at a 
large constant flight-path angle (y(t) = y Q , force of gravity being 
neglected) , the maximum deceleration obtained is directly proportional 
to the velocity and rate of descent. This proportionality can be 
shown to be more general than the limitations imply. For all trajec- 
tories, including those initiated at small flight-path angles, the 
condition of maximum deceleration was derived an equation (l8) . When 
V sin 7 is replaced with h, equation (l8) can be arranged in the form: 


^ an ^max 


—-ft (1 + m 

2gH sin a\ V^/ 


( ) 

v Xi 'max 


-hV 

2gH 
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The approximate expression is comparable to that given (with some 
^^ , ^' er ' ence symbols) in reference 1 and the limitations may be seen 
to be that a = 90° since for high-drag low-lift vehicles, in general, 


2gH 

V2 


<0.01 when the maximum deceleration occurs. 


A special case of interest is the trajectory where the pilot 
maintains a constant rate of descent during the entry. It is presumed 
that the vehicle is initially at an angle of attack of 90° at 
350,000 feet and has a velocity of 25,863 feet per second. The condi- 
tion of constant rate of descent is imposed. Since velocity is speci- 
fied, the initial flight-path angle determines an initial magnitude of 
h and the angle of attack is then varied such as to maintain the 
condition h(t) = h Q as long as it is aerodynamically practical. 

The development of the equations used to calculate this type of 

entry is given in appendix B.^ By using these equations, time histories 
of constant-rate-of-descent trajectories axe presented in figure 18 for 
four different rates of descent determined by the four intial flight- 
path angles selected. 

The data indicate that throughout most of the trajectory only 
moderate changes in angle of attack are required for this type of entry. 
It is further shown that maximum deceleration occurs at the same velocity 
(~13>100 ft/sec) and about the same altitude (*=178, 000 feet) regardless 
of the rate of descent chosen and that therefore the magnitude of the 
maximum deceleration is a linear function of rate of descent (or flight- 
path angle) only. 

A summary plot of the variation of range and deceleration with 
entry angle is shown in figure 19 for the constant h trajectories . 

The linear variation of deceleration with initial flight-path angle is 
given by 


( &n )max 

where y Q is in degrees. Range follows a power relationship to initial 
flight-path angle given by 


^By the use of certain approximations, analytical solutions may be 
obtained which are valid for a large part of each trajectory. The 
approximations and the closed-form solutions are given in appendix C. 
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As a 7 0 _< ^*®7 


but the significance of this fact is not evident. 


Angle of Attack Controlled as a Function of Temperature 

Normally, in order to control the deceleration, it would appear to 
be advantageous to observe deceleration and vary the angle of attack 
accordingly. It has been shown, however, that because of the time 
required to change the flight path a certain amount of anticipation 
should be employed in the control of angle of attack and that a n , 7 , 

or h Q can be used to give this anticipation. A basic trajectory 

variable which may also contain the information useful to the control 
of both heating and deceleration is the skin temperature. This 
variable, unlike a n , can be measured directly and in some cases may 

tend to lead the deceleration and provide the required anticipation. 

Some of the quantities pertaining to the skin temperature obtained on a 
vehicle during an entry are illustrated in figure 20. 

Shown in figure 20 are the calculated data for a typical time 
history of an atmospheric entry. The altitude at the start of the t ime 
history is 550,000 feet and the conditions at 150 seconds are 
h = 550)000 ft, V = 25,863 ft/sec, 7 = - 3 °) and a = 90 °. The angle 
of attack was held fixed at 90 ° during the entire trajectory. The 
quantities shown are: deceleration, a n j time rate of change of 

deceleration, a^ skin temperature for a unit area of material having 
a heat-capacity factor of 0 . 8 , Tp, and material having a heat-capacity 
factor of 0 . 2 , Tp* j the free-stream temperature, Tj the rate of change 
of the two skin temperatures, Tp and Tp*j and the heating rate to the 
skin from the boundary layer, (Q/A) aero . 

It may be seen that Tp and have, in general, the same shape 

and each reaches a maximum at about the same time. Likewise, a n and 
Tp have the same shape , both go negative at the same time although Tp 
leads by about 20 seconds in reaching a maximum. The aerodynamic 

heating rate appears to be strongly influenced by the time rate of 
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change of density or dynamic pressure^ since it also reaches a maximum 
with a^. 


When the insulated unit area of surface material is changed by using 
a different material or a different thickness such that its heat capacity 
is reduced by a factor of 4 (WpC/A reduced from 0.8 to 0.2), the tem- 
perature Tp* increases more rapidly and reaches a higher maximum value 

sooner than the material with the higher heat capacity. Therefore, by 
adjusting the heat capacity, the temperature can be made to lead the 
deceleration and to reach a maximum at about the same time as a^. By 

using this lead, it was thought that Tp* might exhibit characteristics 
similar to a n in the control of angle of attack. Trajectories covering 
the range -of -entry angles were calculated by using AT-,* (base value was 
Tp* at 550,000 feet) to control angle*of attack and the results are 

shown in figure 21. It was found that the trajectories were almost iden- 
tical to those obtained with the a « a n control of figure 7 . At an 

entry angle of 5° the lead effect held the deceleration to 4g but large 
oscillations in deceleration still occurred as they did with the a n con- 
trol. A further disadvantage in the use of Tp or Tp* lies in the 

practical uncertainty of choosing a gain for a variable which does not 
start from a firmly established base value. On the other hand, the fact 
that T^ has about the same variation as ^ during an entry makes Tp 

appear promising as a variable for controlling angle of attack. 

One further combination of variables which appears to offer some 
promise is the quantity Tp* - Tp. As can be seen from figure 20, the 

difference Tp* - Tp is zero at the start of the entry, builds up to a 
positive maximum before and becomes negative at about the same 

time as a n . A study of a number of trajectories shows this to be a 
general characteristic and that the magnitude of Tp* - Tp decreases 
with decreasing entry angles. 

b Dynamic pressure has the same time history as a n but has the mag- 
nitude q = ^a n where = II lb/sq ft for all trajectories of this 

report. Likewise, q = By neglecting the variation of T, it can 

be shown that (Q/ A ) aer o a l 1 ‘55p-l*05. Thus, the effect of p on the 

heating rate is that of an exponential decrease in amplitude which shifts 
the peak magnitude of the heat-transfer rate ahead of the peak decelera- 
tion in time. 
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CONCLUDING REMARKS 


Numerical computations have been made for the trajectories of high 
drag vehicles entering the earth's atmosphere at large angles of attack 
and pertinent data from these computations have been presented. When 
the trajectories were calculated it was assumed that angle of attack 
could be perfectly controlled as a function of some measured quantity, 
held constant at some specified value, or moved (in a continuous manner) 
to maintain one of the other variables at some constant value. The 
results represent large-scale motions of a vehicle which is capable of 
producing lift while descending into the earth's atmosphere from 
initially small flight-path angles and near-satellite velocities. 
Newtonian flow was assumed in calculating the aerodynamic force on the 
vehicle through the trajectory. 

The velocity, altitude, flight-path angle, and angle -of -attack 
requirements necessary to maintain (but not to establish) constant 
deceleration entries of lg, 2g, Jgj and 4g were calculated and used as a 
guide for determining methods of controlling the angle of attack so as 
to limit the deceleration to the arbitrarily chosen value of 3g* When 
the angle of attack was changed from 90 ° in proportion to the decelera- 
tion only, it was found that deceleration could be limited to the 
vicinity of 3g as long as the initial flight-path angle was under 2° 
but at 3° the control was inadequate and large oscillations occurred 
about this mean value of deceleration. 

When the angle of attack was controlled as a function of and 

a n from its initial value of 90 °, deceleration could be held to a 

maximum of 3g for all entry angles to 3°. Trajectories where angle of 
attack was a function of a n only likewise had about the same maximum 
value of deceleration at all entry angles but the level of this maximum 
deceleration was larger (about 4.7g) than when angle of attack was a 
function of both and a n . 

Range and deceleration data were presented for trajectories where 
angle of attack was held constant at values of 84°, 79°, and 74° and for 
trajectories where angle of attack was changed from these values in 
proportion to a n . It was found that, when the entry was initiated at 
these lower angles of attack and then varied as a function of a n , 
smaller and less rapid changes in angle of attack were necessary to 
achieve the desired levels of deceleration than were necessary when a Q 

was 90°. 

The angle -of -attack requirements necessary to maintain constant 
rates of descent throughout the entry were investigated and were found 



27 


to be very moderate. In addition, the maximum deceleration is a linear 
function of that rate of descent only, and the altitude and velocity at 
which the maximum deceleration occurs are virtually independent of the 
rate of descent chosen. Thus, for a vehicle of given wing loading, the 
altitude and velocity at which the deceleration reaches a maximum will 
always be the same regardless of rate of descent. 

By using a simple model of the exposed skin surface of the vehicle, 
the heating rates and skin temperatures were calculated for a number of 
these trajectories. It was shown that, when the angle of attack is 
reduced from 90 ° to 79° or is varied as a function of while the 

deceleration is controlled, the maximum rate of aerodynamic heat 
transfer to the skin is reduced, the maximum skin temperature of the 
vehicle is virtually unaffected, and the total heat to the vehicle is 
slightly increased. The increases in total heat transfer can be minimized 
by holding the maximum desired deceleration for as much of this trajec- 
tory as possible since total heat increases with the total time to enter. 

The feasibility of using a small insulated section of the skin for 
temperature measurements and the use of the measured temperatures for 
angle -of -attack control was investigated. Skin temperature had the same 
time -history characteristics as deceleration, and the time rates of 
change of these two variables were likewise very similar throughout the 
trajectory. Temperatures could be made to lead the deceleration by 
reducing the heat capacity of the material. With and without this lead, 
the use of temperature to control angle of attack produced trajectories 
very similar to those obtained when angle of attack was controlled by 
deceleration. A simple means of anticipating the increase in decelera- 
tion and heating rates is to measure the difference in temperature 
between two adjacent sections of the skin of the vehicle having different 
levels of conductivity. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., October 13, 1958* 
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APPENDIX A 

EQUATIONS FOR CONSTANT DECELERATION TRAJECTORIES 


The basic differential equations for an entry into the earth's 
atmosphere are 


V 

— = -a n sin a - sm 7 
& 


(Al) 


For 



a n cos cx "t* 



h = V sin 7 

constant-deceleration trajectories, it 


a n 


-a z 


= e - h / H V 2 

2W7S V 


cos 7 


(A 2) 
(A5) 


is further required that 


= Constant 


(A4) 


The object is to find the variation of angle of attack which will satisfy 
these equations at every point on the trajectory. The other variables 
are dependent on the proper choice of a and they are determined in the 
process . 


From equation (a 4) the relationship between altitude and velocity 
for the constant values of deceleration and wing loading is given by 


h 


- H l0 s e 


°FPe v£ 
2W/S a n 


(A5) 


The condition that a n be zero (see eq. (17)) defines a necessary rela- 
tionship between h and V; namely, 


Vh = 2HV 


A second differentiation gives 


h 



(a6) 


(A?) 


Thus h and its time derivatives are defined in terms of V and its 
derivatives . 
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The flight-path angle may likewise be related to velocity. From 
equations ( A 3 ) and (a 6 ) , 


7 


sm 


-1 


h 

V 


7 


sm 


-1 2H 


V 


Differentiating equation (A 8 ) and making the substitution 


yield 


cos 7 


1 




1/2 


7 


2H 

V cos 




2 


(A 8 ) 


(A9) 


(A10) 


The time rate of change of velocity is derived from equations (Al) 
and (A 8 ) in terms of angle of attack and velocity as 


-ga n sin a 



(All) 


Differentiation of equation (All) with respect to time and elimination 
of V give an expression for the second derivative of velocity 


V 


. 2 
sm 7 


Ha n cos an/ 


1 + 


2gH 


(A12) 


Equations (A2), (A 8 ) , (A10) , (All), and (A12) together with an equation 
for a represent six equations with six unknowns at each value of time 
The equation for a is unspecified but one may be manufactured in sev- 
eral ways, the most convenient being 


a 



(M3) 


Thus for a given change in velocity, equations (All) and (A13) depend 
only on the corresponding change in angle of attack. The integration 
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process enters the equations only through calculation of the dc term, 
a variable which is generally small. Since a constant-deceleration 
entry is being considered, V should be very nearly constant over each 
A interval of integration. In the calculations of this paper, for each 
small reduction in velocity from V n to V^, equation (A13) was 

applied in the form 

an " v n - v n-l 2(57.3) \ n + Vn_1 ) (M4) 


to determine the nth value of a. The corresponding value of time is 
given by 


^n ~ t n-l + 


n 


'n -1 


§( v n + V,) 


(A15) 


It is useful in determining the initial velocity to combine equa- 
tions (A2) , (A 8 ) , (AlO) , (All) , and (A12) into the following equations : 


2 sin q tan 7 



V^ j cos 7 
g r / a n 


cos a 


V tan y • 

x a 

a n g tan a 


(ai6) 


sin 7 


^gHaj^ sin a 
V 2 + 2gH 


(A17) 


Thus, when a n , a, and d are specified, these two equations determine 
the corresponding value of the velocity. 

For a given deceleration the maximum entry angle is obtained when 
a, is 90°, as may be seen from equation (A17) . Since the terms involving 
a in equations (A12) and (Al 6 ) become zero when a is 90°, this condi- 
tion affords an initial condition where all the variables are known. 
Thereafter, a AV interval is- chosen, a value of a is assumed at the 
new velocity, and equations (All), (Al4) , (A12) , (A 8 ) , (AlO) , and (A2) 
are solved in that order to yield a new value of a. An iteration 
process is used until the assumed and derived values of a are in 
agreement. The iteration process is easily handled on a high-speed 
digital computer and at the higher values of velocity the process is 
highly convergent. 


A -very good approximate closed-form solution is available based on 
the assumption that 
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sin a = Constant =1.0 

With this assumption, equation (All) may be integrated to give the veloc- 
ity as a function of time 

V(t) = Jp(t ) 2 + 2gH - 3 (t) (A18) 

where 

P(t) = i a n gt - V G + 

^ v o 


When V is given as a function of the time and sin a is treated as a 
constant with a value of unity, the other variables may be directly cal- 
culated as follows: 
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(A19) 


Making the indicated substitutions in equation (A2) or applying the 
assumption to equation (Al 6 ) then yields an expression for a(t) 


cos a 


2 tan 7 



VMcos 7 
SrJ a n 


(A20) 
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This approximate solution leads to results very close to those obtained 
by the iterative integral method as long as a remains near 90°. How- 
ever , the solution becomes imaginary when the quantity 

2gHa n 

“P > 1 

+ 2gH 

At this velocity no real value of y exists which will satisfy the 
equations . 
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APPENDIX B 

EQUATIONS FOR CONSTANT RATE OF DESCENT TRAJECTORIES 


The basic equations for trajectories into the earth’s atmosphere 
where rate of descent is maintained at some constant value are given by 



V 

— = -a n sm a - sm 7 

0 

(Bl) 

V 7 . 

g 

/v 2 \ 

' = a n cos a + — - 1 cos 7 

n \gr j 

(B2) 


h = h 0 = V sin 7 

(B3) 


, Ve -h/H 2 

n 2 W/S 

(B4) 

From equation (B 3 ), the 

following variables are extracted: 
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(B 5 ) 


h = h 0 t + h Q 

(b6 ) 


. -1 ^0 
7 = sm — 

(B 7 ) 


• -V tan 7 

7= V 

(b8) 

Equation (b 8) expresses the required relationship between 7 and V 
necessary to maintain a constant rate of descent and relates equation (Bl) 
to equation (B 2 ) . By substituting equation (B8) into equation (B 2 ) to 
eliminate y and then using equation (Bl) to eliminate ^ , the following 
relationship between 7 and a is obtained: 

a n 

2 

cos(a +7) = 1 - — cos^7 
gr 

(B9) 
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Therefore if V and 7 are to be related properly, then a must fol- 
low the relationship 


a 


-7 + cos 


-1 


1 A v 2 2 \ 

— 1 cos 7] 

a n \ gr 7 


(BIO) 


If the initial values of h, V, and 7 (or h) , are specified, 
the other variables may be calculated at these initial conditions. An 
integration process is used to establish approximate values of V and 
7 at the end of the given time interval. Equations (b6), (B 4 ), (BIO), 
(Bl), and (B8) are solved in that order and a new approximation to V 
and 7 made. The process may be repeated but convergence is generally 
rapid. From the calculated values of V and 7, equation (B3) is solved 
for h at that point on the trajectory. Any inaccuracy in the integra- 
tion process will produce an h different from its initial value and 
this error is carried over to the next point on the trajectory. Thus 
the values of h calculated at each point along the trajectory depend 
directly upon the cumulative accuracy of the integration process and 
serve as a measure of that accuracy. The effect of these errors shows 
up in. figure 18 as a slight departure of the rate of descent from its 
initial value. In all cases this error was less than 3.5 percent at 
the end of the calculations. 
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APPENDIX C 


APPROXIMATE ANALYTICAL SOLUTIONS TO CONSTANT -DECELERATION 
AND CONSTANT -RATE -OF -DESCENT TRAJECTORIES 


It has been pointed out that, by using the ?! Z M function derived 
by Chapman (ref. 2), certain closed-form solutions can be made for the 
trajectories of constant deceleration and constant rate of descent. 
Solutions can also be obtained directly from the basic equations of this 
paper by applying the approximations used in reference 2 , and this method 
is used here. In so doing, it is felt that a better physical under- 
standing of the origin and limitations of each solution can be retained. 


the 


For most of these closed-form solutions, it is 
following assumptions: 

(a) [7 | is small; thus, cos 7=1 and tan 

is small; thus, sin a = 1 and 


(b) 

(c) 


* - a 

2 


2 gH 


»> 1 


necessary to make 

7 = sin 7 
cot a = cos a 


(d) 



a 


n 


Although these assumptions may appear to be rather restrictive, they are 
generally valid on that portion of the trajectory of the most interest, 
namely, the region of maximum heating and deceleration. It may be seen 
from equation (Al) that assumption (a) is violated at the beginning of 
the entry when a n ~ -7 but the effect is not serious since V/g is 
small at this point. At the later stages of the entry all the assump- 
tions are generally invalid. 

Chapman’s symbol for the ratio of velocity along the trajectory to 
the initial circular satellite velocity 


u 


V 

\fer 


is used in this appendix. Thus, for entries from circular and near- 
circular orbits, u Q can be set to unity. 
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Constant -Deceleration Trajectories 

Closed-form solutions for constant -deceleration trajectories based 
only on the assumption that sin a = 1 are given in the latter part of 
appendix A. Further simplifications, based on the foregoing assumptions, 
lead to the analytical solutions which follow. 


a n , its compo- 


In view of the assumptions ^ the normal acceleration _ 

nent in the direction of the velocity a n sin a, and its horizontal com- 
ponent a n sin a cos 7 are all approximately the same. For a specified 

constant value of a n , equation (A19) gives the expression for the flight- 
path angle 


sin 7 


n 




u 4501? 


From equation (A20) the angle of attack becomes 


(Cl) 


t - 2 a ... 

f'nfi ri m ± & n ( 

1 - u 2 



D r — 2 

— u 
2H 

a n 


(C2) 

_ 1 - u 2 

a n 



a n 

225 U 2 



in time from u 

n -1 

to u n 

is 

calculated from 

At . Mfe , 

S a n a „ l n -1 

- a n) 


(C3) 


The incremental change in range is given by 



1 

t n 

r n 

— 


/ V dt 

vr / 

r 0 

6 n _l 



r Vn 


- X 

/ V dV 


rV 

Vi 



V U n-l “ u n) 


(C4) 
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The altitude-velocity relation given by equation (A4) becomes, in terms 
of u: 


h 


H log e 


C F p e 

2 W/S 


gr 



(C5) 


When desired, the derivatives of the various variables are given 
with fewer restrictions by equations (A19). 


Constant -Rate -of -Descent Trajectories 

From equation (B 7 ) the flight-path angle during the descent is 
given by 

V sin 7 
sin 7 = — M 

V 


or 


u o . 

sm 7 = — sin 7 C 
u 


(c6) 


The deceleration can be obtained in closed form by integrating 

h = V sin 7 = V Q sin y Q (C7) 


Since 


and 


h =^v 
dV 


-V 

g 


equation (C7) becomes 


-ga n dh = V 0 sin y Q dV 

Dividing both sides by V 2 , substituting for a n its definition given 

by equation (B4) , and integrating from initial to some arbitrary condi- 
tion of V, h gives 
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-g 


/ C FPe\ 

1 2W/S y j 


-h/H 

e dh = V Q sin 7 ( 


r v dv 

V 2 


gH 




"(Vo Sin 7o)v[l - £ 


and finally in terms of u 


a n - H r ( u c 



u 

)u 1 - — 

+ — 

V l uj 

' u^ 


o 

sin 7 dull 

-P 

V 

u ol 


ln o 


(c8) 


This deceleration is a maximum when 


u r 


u = 


(C9) 


and the value -of this maximum deceleration is 

Nmax ~ - 22 5 a o(*o sin 7 0 ) 




(CIO) 


where j Q is small and is given in degrees. 

Since a n is defined, the variation in angle of attack is given by 
writing equation (B9) in "the form 

1 ( p p N 

cos a cos 7 - sin a sin y = — (1 - u cos y) 

a n 


By virtue of the assumption that cos y = sin a = 1, this equation can 
be approximated with 


— ~ sin 7 4- — (l - n 2 ) 


cos a ~ ~ sm y 

D 


(Cll) 
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or as a function of u only 


cos a - h « U ° Sin y o 
D u 


900 (u Q sin 7 0 )u 


(C12) 


The altitude is given by equation (B6) in terms of known ini- 

tial condition h Q , h Q and time* Time is related to u by integrating 

the approximate expression 


- ~ -a r 


900 ( 


(u D sin 7 o)u 

(1 -J-) 

\ / 

\ u o/ 


Thus 


J n 


dt 


^n-1 


up r Un du 

900g(u o sin 7o ) u(u Q - u 


yields 


At = fc n - t n-l * r 


0.9 i ^^o ~ u n-l) 

— l0f 5e r-f 

u o sin 7 0 u n-l( u o - u n) 


The incremental change in range is given by 


A ( 


* i r tn v dt - r Un sjs 

^ rJ t n . x J u n _ L - a n 


(C13) 


Substituting equation (C8) for a„ and integrating yield 



1 

900 sin 7 q 


log e 


u o - u n-l 
^o - ^n 


(C14) 


Some of the deductions and comments of the text are more clearly 
indicated by the foregoing approximations. For example, the deductions 
that, for a constant rate of descent, the velocity at maximum decelera- 
tion is independent of entry angle and that the solutions are independ- 
ent of W /CpS and density, except for the altitude -velocity relation- 
ship, are seen in these solutions* 


■# 
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All the foregoing equations may be derived in the nondimensional 
form of reference 2.^ In the notation of reference 2, the horizontal 
component of deceleration is given therein as 

a 0 » P rlIZ ~ 30uZ 
cos cp 


where B = — of this paper, u = —S— = — - S . 9^ 2 . is the nondimensional 
H |gf Jgr 

horizontal component of velocity , cp = 7 of this paper } and the function 
Z (in terms of the notation of this paper) is 

Z = JgH e”* 1 /^ V sin a cos 7 

2W/S V 



sin a cos 


7 


For sin a = cos 7=1, 


30uZ ~ a n 

which is the reason the results for the constant-deceleration trajectories 
are the same in either system of coordinates as long as the assumptions 
are valid. 


J-For constant-deceleration trajectories Z = and for constant- 

30u 


rate -of -descent trajectories Z = -30u o sin cp 


- €)• 
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Figure 3-- Trajectory data for entry into the earth's atmosphere at a 
constant value of deceleration. Time histories are shown for decel- 
erations of lg, 2g, 3g, and 4g. In each case, a 0 = 90°. 
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Figure 7.- Comparison of trajectories controlled with a = 90° - 
when initial flight-path angles are - 1 , -1°, -2°, and -3°. 

V Q = 25,863 feet per second; ho = 350,000 feet. 
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(a) 7 0 = - 1 °. (b) 7 0 = - 2 °. (c) y Q = - 3 °. 

Figure 12 .- Comparison of trajectory data, skin temperatures, and aerody- 
namic heat-transfer rates for initial flight-path angles of -1°, -2°, 
and - 3 ° when trajectory is controlled with a = 90 ° = Constant and 
a = 90° = 250^. V D = 25,863 feet per second; h Q = 350,000 feet. 
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Figure 18.- Trajectory data for entry into the earth's atmosphere at a 
constant rate of descent. Rates of descent determined from initial 
flight-path angles: V n = 26,000 feet per second; h Q = 350,000 feet. 









Time, sec 

Figure 20.- Comparison of deceleration temperature and aerodynamic- 
heating data calculated for a typical atmospheric entry, 
ho = 551,513-5 feet; V Q = 25,633-5 feet per second; y D = -2.96°. 

(At 150 seconds, h = 350,000 feet; V = 25,863 feet per second; and 
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